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FILTERING FOR SPATIAL AUDIO RENDERING 
Background 

Accurate spatial reproduction of sound has been 
demonstrated to significantly enhance the visualization of 
5 three-dimensional (3-D) multimedia information, 

particularly with respect to applications in which it is 
important to achieve sound localization relative to visual 
images. Such applications include, without limitation, 
immersive telepresence; augmented and virtual reality for 
10 manufacturing and entertainment; air traffic control, pilot 
warning, and guidance systems; displays for the visually/or 
aurally- impaired; home entertainment; and distance 
learning. 

Sound perception is known to be based on a 
15 multiplicity of cues that include frequency- dependent level 
and time differences, and direction-dependent frequency 
response effects caused by sound reflection in the outer 
ear, cumulatively referred to as the head-related transfer 
function (HRTF) . The outer ear may be effectively modeled 
20 as a linear time -invariant system that is fully 

characterized in the frequency domain by the HRTF. 

Using immersive audio techniques, it is possible to 
render virtual sound sources in 3-D space using an audio 
display system, such as a set of loudspeakers or 
25 headphones. The goal of such systems is to reproduce a 
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sound pressure level at the listener's eardrums that is 
equivalent to the sound pressure that would be present if 
an actual sound source were placed in the location of the 
virtual sound source. In order to achieve this result, the 
5 key characteristics of human sound localization that are 
based on the spectral information introduced by the HRTF 
must be considered. The spectral information provided by 
the HRTF can be used to implement a set of filters that 
alter nondirect ional (monaural) sound in the same way as 

10 the real HRTF. Early attempts at the implementation of 

HRTFs by filtration were based on analytic calculation of 
the attenuation and delay caused to the soundfield by the 
head, assuming a simplified spherical model of the head. 
More recent approaches are based on the measurement of 

15 individual or averaged HRTF' s that correspond to each 
desired virtual sound source direction. 

In addition to simulating the effects of cues that 
operate on the human ear, effective spatial audio rendering 
engines must also accurately simulate the virtual ambient 

20 in which the listener is to experience the spatially 

reproduced sound. To this end, a spatial audio rendering 
engine typically retrieves a set of reverberation paths 
that extends between the sound source and the listener. 
Reverberation paths may be retrieved in accordance with a 

2 5 number of known techniques, prominently including beam 

tracing. Using the reverberation paths, the spatial audio 
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rendering engine then synthesizes a signal that faithfully 
replicates an actual listening experience. 

Heretofore, realization of the above-described process 
in a manner that results in a convincing audio simulation 
5 has been found to be computationally daunting. 

Accordingly, what is required is an approach to spatial 
audio rendering that, in one regard, reduces computational 
complexity, while concurrently affording the desired degree 
of simulation quality. In another regard, there exists a 
10 need to provide an audio rendering engine that admits of a 
capability to effect a counterbalance, at a user's 
discretion, between computational complexity and quality of 
audio reproduction . 

Brief Description of the Drawings 
15 The subject spatial audio rendering technique may be 

better understood by, and its many features, advantages and 
capabilities made apparent to, those skilled in the art 
with reference to the Drawings that are briefly described 
immediately below and attached hereto, in the several 
20 Figures of which identical reference numerals (if any) 
refer to identical or similar elements, and wherein: 

FIG. 1 is a graphical representation of the manner in 
which the physical characteristics of a virtual audio scene 
may, in one embodiment of the invention, be considered in 
25 the design of a spatial audio rendering system. 
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FIG. 2 is a block diagram of a spatial audio rendering 
system in accordance with an embodiment of the invention. 

FIG. 3 is a block diagram of an exemplary processor- 
based system into which embodiments of the invention may be 
5 incorporated . 

Skilled artisans appreciate that elements in Drawings 
are illustrated for simplicity and clarity and have not 
(unless so stated in the Description) necessarily been 
drawn to scale. For example, the dimensions of some 
10 elements in the Drawings may be exaggerated relative to 
other elements to promote and improve understanding of 
embodiments of the invention . 



Detailed Description 
15 Referring now to FIG. 1, depicted therein is a 

generalized representation of the methodology according to 
which, in at least one embodiment of the invention, the 
physical characteristics of a virtual audio scene may be 
captured and quantified so that spatial audio rendering may 
20 be effectively implemented. The intended result of a 
spatial audio rendering system is to reproduce (or 
simulate) the listening response of a human being at a 
defined position in a virtual scene. The virtual scene be 
presented, for example, in the application of computer 
2 5 graphics, music playback, sound tracks, and other 

entertainment content. The listening experience is 
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understood to be a function of the sound sources and the 
ambient scene geometry and material properties. 

Essentially, the spatial audio rendering system 
operates to capture each reverberation path that couples a 
5 sound source to listener. Generally speaking, a 

"reverberation path" may be here understood to be a trace 
that represents sound propagation in a scene by taking into 
account interaction with a single or multiple obstacles. 
In this regard, beam tracing may be used as a technique for 

10 modeling the interaction of sound with obstacles. In 
general, the beam tracing approach assumes specular 
reflection of sound beams off relevant obstacles. Simple 
geometrical calculations allow the definition of 
reverberation paths from the sound source to the receiver 

15 point. Reverberation paths may be represented 

geometrically as a form of polyline. Source image 
positions are calculated that constitute each real source 
in the scene. As result, the scene, which contains a 
number of obstacles and real sources, is represented as 

2 0 free space that contains a set of source images and a 
receiver. Equivalently polyline beams emitted by real 
source and received by receiver are replaced with set of 
source images, each source image emitting one linear beam 
(a reverberation path) received by a receiver. 

2 5 A reverberation path may be said to be "captured" by 

virtue of mathematical characterization in terms of, for 
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example, the attenuation and delay imparted to a signal 
source by the reverberation path. Accordingly, for signal 
processing purposes, each reverberation path may be 
represented by a filter that imposes a predetermined 
5 frequency-dependent attenuation on a source image signal. 

Filters corresponding to the respective reverberation 
paths are coupled to the signal source (s) to generate a 
reverberant signal that is associated with each 
reverberation path. The reverberant signals are then 

10 accumulated to produce a resultant (simulated) signal. The 
resultant signal is delivered to the listener through an 
audio display system, e.g., headphones or loudspeakers. 

As graphically represented in FIG. 1, each 
reverberation path may be traced and represented as an 

15 source image that is characterized, according to the 

geometry of the virtual scene, by a set of coordinates, 
e.g., azimuth and elevation. As indicated above, a source 
image technique is used to model sound propagation and 
interaction with obstacles. Once specular sound 

2 0 reflections are assumed, a scene containing obstacles may 
be simulated by free space containing real and 
corresponding source images. 

Consider, for example, a scene with one reflecting 
wall. In this case there exists a direct sound propagation 

2 5 path from source to receiver, as well as one reverberation 
path from the source to the wall and from the wall to the 
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receiver. The wall may be considered in the nature of 
mirror. Therefore, the real scene (containing a wall) may 
be simulated by free space containing a real source and a 
mirrored (image) source. The foregoing constitutes the 
5 essence of the source images construct, as applied to 
spatial audio rendering. 

Be aware, however, that the above example illustrates 
a first-order source image that models sound interaction 
with a single obstacle. In general, a scene may contain a 

10 greater number of obstacles, and the order of reflections 
(source images) is then concomitantly much higher. A 
second order source image can be calculated by mirroring 
the first order source image in another obstacle. That is, 
a second-order source image models sound propagation from a 

15 source to a receiver, and includes interactions 
(reflections) with two obstacles, etc. 

Material properties, such as frequency- dependent 
reflection coefficients, of the virtual scene are also 
relevant and are considered in the design of filters 

2 0 employed to characterize a given respective reverberation 

path. In this form, the characterization process enables a 
specific filter design, specified by filter coefficients, 
that corresponds to each reverberation path. (The manner 
in which the filters are designed is not considered here to 

2 5 be an aspect of the present invention. Suffice it to say 
that practitioners skilled in the art of digital signal 
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processing techniques possess expertise adequate to 
synthesize digital filters that implement frequency- 
dependent amplitude and delay characteristics. See, for 
example, D. Schlichtharle , "Digital Filters: Basics and 
5 Design," Springer, 2000.) 

As represented in FIG. 1, the characterization process 
results in, for example, a set of filter coefficients that 
correspond to each reverberation path. A filtering module, 
designed in accordance with the coefficients, accepts an 

10 input signal that originates with a sound source and 

filters the signal according to the parameters of the set 
of source images that correspond to reverberation paths. 
As indicated above, in one embodiment, filtering comprises 
the application of a frequency- dependent attenuation factor 

15 and the insertion of a time delay. The reverberant signals 
(filtered source image signals) are accumulated to 
synthesize a resultant signal. Typically, the resultant 
signal is divided into at least two channels, e.g., left 
and right; although in alternative embodiments, more than 

2 0 two channels may be created. The output channels may then 
be applied to one or more audio display systems, such as, 
for example, a loudspeaker system or a headphone system. 

In alternative embodiments, prior to the accumulation 
of the reverberant signals and delivery of an output signal 

25 to the audio display system, an additional filter (which 
may be considered a "post- filter" ) may be applied. The 



characteristics of the post-filters are dependent on the 
nature of the audio display system and are also dependent 
on the coordinates of the source images. For example, as 
indicated above, HRTFs may be applied to the reverberant 
5 signals prior to accumulation and application to a 
headphone system. 

In addition, in applications where a loudspeaker 
system is incorporated as an audio display device, 
filtering appropriate to the Ambisonic technique may be 

10 applied. As is known to those skilled in the art, in the 
application of the Ambisonic technique, an output signal 
for each speaker is produced as a weighted sum of 
individual reverberation path signals. The weight 
coefficients may be calculated from the source image 

15 coordinates and loudspeaker layout. 

Ambisonic sound processing is a set of techniques for 
recording, studio processing and reproduction of the 
complete sound field experienced during the original 
performance. Ambisonic technology decomposes the 

2 0 directionality of the sound field into spherical harmonic 
components. The approach uses all speakers in a system to 
cooperatively recreate these directional components. That 
is to say, speakers to the rear of the listener help 
localize sounds in front of the listener, and vice versa. 

25 Ambisonic decoder design aims to satisfy simultaneously and 
consistently as many as possible of the mechanisms used by 
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the ear/brain to localize sounds. The theory takes account 
of non-central as well as central listening positions. In 
an Ambisonic decoder, the spherical harmonic direction 
signals are passed through a set of shelf filters that have 
5 different gains at low and high frequencies, wherein the 
filter gains are designed to match the panoply of 
mechanisms in which the ear and brain localize sounds. 
Localisation mechanisms operate below and above about 700 
Hertz (Hz) . The speaker feeds are then derived by passing 
10 the outputs from the shelf filters through a simple 

amplitude matrix. A characteristic of Ambisonic decoder 
technology is that it is only at this final stage of 
processing that the number and layout of speakers is 
considered . 

15 For a thorough understanding of the subject spatial 

audio rendering technique, refer now to FIG. 2, which is a 
block diagram of a spatial audio rendering system 20 that 
is implemented in accordance with one embodiment of the 
invention. As illustrated in FIG. 2, system 20 comprises 

2 0 an input stage 211 that may be coupled to an audio input 

signal source 210. In one embodiment, audio input signals 
that are stored in, or are transmitted from, signal source 
210 may be provided as digital files, such as, for example, 
AFFI, WAV or MP3 files. However, the scope of the 

25 invention is not constrained by the nature of input files, 
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and embodiments of the invention extend to all manner of 
digital audio files, now known or hereafter developed. 

Input stage 211, in one embodiment, may be constructed 
to divide the digital audio input signal into a number of 
5 timewise-overlapping windows. 

There exist numerous techniques to divide a time- 
domain input signal into windows. The primary purpose of 
the signal windowing is a further calculation of the 
frequency domain signal spectrum, which may be 

10 accomplished, for example, using a Fast Fourier Transform 

(FFT) . 50% overlapped sinusoidal windows may be typical in 
one embodiment of the invention. The length of the window, 
in one embodiment, may vary from 256 to 2 048 samples of the 
input time-domain signal. Other arrangements of the 

15 window, including the overlapping ratio and length are also 
possible. Skilled practitioners, in the judicious exercise 
of a designer's discretion, may select window shape, 
overlapping ratio and length to obtain more nearly optimal 
results that are tailored for an individual application. 

2 0 However, window shape, overlapping ratio and window length 
are not constraints on the scope of the invention. 

The output of input stage 211 is coupled to an FFT 
(Fast Fourier Transform) module 212. In a manner well 
understood by practitioners acquainted with digital signal 

25 processing (DSP) techniques, FFT module 212 operates to 

transform each of the timewise-overlapping windows created 
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by input stage 211 to a frequency- domain equivalent, that 
is, into a frequency- transformed window. The frequency- 
transformed windows are stored in a cyclic input buffer 
214. In practice, cyclic buffer 214 comprises a number of 
5 distinct buffers 214a, 214b, 214n, each of which stores 
one of the frequency- transformed windows. The length of 
the buffers may be designed to correspond to the length of 
the longest delay interposed by a reverberation path. In 
general, a buffer adequate to insert a delay of one second 

10 (at the applicable system clock rate) may generally be 

sufficient, although other implementations would suggest 
different buffer sizes. 

A spatial audio rendering engine 216 may be 
constituted from a plurality of source image processing 

15 kernels 216a, 216n . In the manner indicated in FIG. 2, 
each of the source image processing kernels may be 
selectably coupled (as described below) to an output of one 
of the cyclic input buffers 214a, 214n. Coupling of an 
input buffer to one of the source image processing kernels 

2 0 may be effected under software control, for example. 

In addition, and as depicted in FIG. 2, in operation, 
each of the source image processing kernels 216a, 216n 
is also associated with one of the filters 215a, 215n 
that constitute filter bank 215. Filters 215a, 215n are 

25 constructed, as described above and depicted in FIG. 1, to 
characterize the reverberation paths alluded to above. 



That is, each of the filters 215a, 215n in filter bank 
215 is designed to impart to a source image a frequency- 
dependent attenuation that simulates a reverberation path. 
Accordingly, each of the filters 215a,..., 215n corresponds 
5 to a reverberation path. Filters 215a, 215n may be 

realized as digital filters having characteristics that are 
defined by predetermined filter coefficients. 

In one embodiment, source image processing kernels 
216a, 216n operate in the following manner, under 

10 software control, for example, to process selected ones of 
the frequency- transformed windows stored by cyclic input 
buffer 214. Specifically, in one embodiment, for each 
reverberation path that has been identified with respect to 
a virtual scene, a signal delay is determined for each path 

15 between a source image and the listener. The delay may be 
determined in accordance with any of a number of 
techniques, such as, for example, by the acquisition of 
empirical data or as a result of a mathematical 
calculation, based on, for example, the distance between 

2 0 the source image and the listener. Software simulation may 
also be employed. Once a signal delay is attributed to each 
reverberation path, the transformed window having a delay 
that is closest to the delay attributed to the 
reverberation path is identified and thereby matched to 

25 reverberation path. In this regard, it should be noted 
that, as a matter to be determined in the judicious 



discretion of the system designer, smaller time-delay 
distances between consecutive frequency- transformed windows 
stored in buffer 214 result in finer granularity in the 
match between reverberation path, i.e., source images and 
5 available transformed windows. However, the improvement in 
matching is acquired at the expense of an increase in the 
number of frequency-transformed windows that must be 
available and, therefore, the number of FFTs that must be 
performed. 

10 In the above -described manner, the transformed windows 

stored in respective ones of the cyclic input buffers 214a, 

214n are selected for concurrent processing by 
associated ones of the source image processing kernels 
216a, 216n. Essentially, the source image processors 

15 operate to apply an appropriate one of the filters 215a, 

215b, 215n to each of the selected transformed windows. 
That is to say, in one embodiment, given a transformed 
window that has been matched to a reverberation path and 
that has been assigned for processing by a source image 

2 0 processing kernel, then processing is performed in 

accordance with parameters established by the filter that 
corresponds to the reverberation path. 

Consequently, the source image processing kernels 
concurrently provide a plurality of output signals, which 

25 may be denominated here as "frequency- domain reverberants . " 
Each of the frequency- domain reverberants corresponds to a 



delayed and attenuated version of a source image that is 
associated with a reverberation path. Delay is effectively 
imparted to a source image operation of the cyclic buffers. 
Frequency-dependent attenuation is imparted by virtue of 
5 the application of a particular filter that has been 

characterized in conformance with the reverberation path. 

In some embodiments, the system may also include a 
table 213 of HRTFs . The table (which may constitute any 
form of suitable storage device) contains a number of HRTFs 

10 that, much like filters 215a, 215n, are matched to an 
source image reverberation path) . Consequently, as 
transformed windows are selectably applied to a respective 
source image processing kernels for processing in 
accordance with appropriately matched filters 215a... 215n, 

15 so too are appropriate ones of HRTFs 213a, 213b, 213n. 
Therefore, in such an embodiment, the reverberant outputs 
of the source image processing kernels represent a delayed 
version of an source image that has been specifically 
attenuated by one of filters 215a...215n to conform to the 

2 0 attenuation interposed by the reverberation path and by 
one of the HRTFs 213a, 213b, 213n to simulate the 
auditory response of a human being to an source image that 
is displayed through headphones. Recall here that HRTFs 
differ as a function of source image coordinates. 

25 Therefore, HRTFs are likewise matched to specific source 
images . 
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As illustrated in FIG. 2, the outputs of the source 
image processing kernels (i.e., reverberants) are, in one 
embodiment, coupled to parallel left (L) and right (R) 
channels 217 and 218, respectively. Each channel comprises 
5 a respective signal combiner (217a, 218a) , output buffer 
(217b, 218b) , Inverse Fast Fourier Transform (IFFT) module 
(217c, 218c) , and interstage buffer (217d, 218d) . 

As to operation, the concurrent reverberant outputs of 
appropriate ones of the source image processing kernels are 

10 coupled to the inputs of the respective left and right 

channel signal combiners 217a and 217b. The outputs of the 
signal combiners, denominated here " frequency-domain 
resultants, " are buffered in respective left and right 
output buffers 217b and 218b and are applied to respective 

15 IFFT modules 217c and 218c. IFFT modules 217c and 218c 

transform the frequency- domain resultant signals into the 
time-domain equivalents, i.e., time-domain resultants. The 
left and right time-domain resultant signals are coupled 
through respective interstage buffers 217d and 218d to an 

2 0 interleave module 219. 

In a manner familiar to those skilled in the art, 
interleave module 219 imparts a standard formatting 
convention that is applicable to the storage and 
transmission of multichannel audio data. For example, with 

25 respect to stereophonic audio data that comprises a Left 

(L) and Right (R) channel, samples are taken in a L, R, L, 



R, L, R, ... sequence. Interleave module 219 operates to 
interleave a sequence of left channel signals (L, L, L, ...) 
and right channel signals (R, R, R, ...) to produce an 
interleaved channel sequence, L, R, L, R, L, R, that can 
5 be stored in a WAV file or played back using a computer 

audio card. The output of interleave module 219 is coupled 
to an audio display device, which may be, for example, a 
loudspeaker system or a headphone set, although other forms 
of audio display devices, now known or hereafter developed, 

10 may be used with the invention. 

The embodiment described immediately above is 
particularly advantageous in applications where the number 
of reverberation paths is relatively small (say, up to 100) 
and relatively fine granularity is required of the source 

15 images. In this context, it is deemed appropriate that the 
input signal, initially provided to the spatial audio 
rendering engine in the time domain, be converted to the 
frequency domain and stored in a cyclic buffer as 
frequency-domain transforms. Consequently, one FFT is 

2 0 required for each channel (e.g. Left and Right) . 

Alternately, the audio input may be coupled directly 
(without FFT) to the cyclic buffer and stored in the time 
domain. Depending on the reverberation path and 
corresponding time delay, the signals stored in respective 

25 buffers are selected and transformed through the 

application of a respective FFT module, so that one FFT 



module is required for each reverberation path. After 
application of the FFT, reverberation path filters, HRTF 
filters and other (if any) filters may be applied to the 
frequency-domain signal. An I FFT is applied to each 
5 channel signal after summation of the individual 
reverberation path signals. 

Furthermore, in some applications the number of 
reverberation paths may be large, greater than 100, for 
example. Specifically, in a small room with complex 

10 geometry and highly absorbent materials, the reverberation 
time is typically quite short, but the number of 
reverberation paths may be significant. Consequently, a 
large number of reverberation paths will share a similar 
delay. In this context, an alternative embodiment may be 

15 warranted in which the use of a matrix filter may be 

invoked. According to the approach, filters corresponding 
to reverberation paths that are matched to the same window 
may be aggregated. As a result, filtration is reduced to 
the multiplication of two matrices of size (M)x(N), where M 

2 0 is the number of windows and N is the length of each 

window. In this embodiment, the computational complexity 
of filtration does not increase with the number of 
reverberation paths. However, when the number of 
reverberation paths is small, the matrix filter is then 

25 only sparsely populated. In this context, the matrix filter 
approach imposes substantial computational overhead. 



FIG. 3 is a block diagram of an exemplary processor- 
based system into which embodiments of the invention may be 
incorporated. With specific reference now to FIG. 3, in 
one embodiment the invention may be incorporated into a 
5 system 300. System 300 is seen to include a processor 310, 
which may include a general -purpose or special -purpose 
processor. Processor 310 may be realized as a 
microprocessor, microcontroller, an application-specific 
integrated circuit (ASIC) , a programmable gate array (PGA) , 
10 and the like. As used herein, the term "computer system" 
may refer to any type of processor-based system, such as a 
mainframe computer, a desktop computer, a server computer, 
a laptop computer, an appliance, a set-top box, or the 
like. 

15 In one embodiment, processor 310 may be coupled over a 

host bus 315 to a memory hub 330, which, in turn, may be 
coupled to a system memory 320 via a memory bus (MEM) 325. 
Memory hub 33 0 may also be coupled over an Advanced 
Graphics Port (AGP) bus 333 to a video controller 335, 

20 which may be coupled to a display 337. AGP bus 333 may 
conform to the Accelerated Graphics Port Interface 
Specification, Revision 2.0, published May 4, 1998, by 
Intel Corporation, Santa Clara, California. 

Memory hub 330 may also be coupled (via a hub link 

25 338) to an input/output (I/O) hub 340 that is coupled to a 
input /output (I/O) expansion bus 342 and to a Peripheral 
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Component Interconnect (PCI) bus 344, as defined by the PCI 
Local Bus Specification, Production Version, Revision 2.1 
dated in June 1995. The I/O expansion bus (I/O EX PAN) 342 
may be coupled to an I/O controller 346 that controls 
5 access to one or more I/O devices. As shown in FIG. 3, 
these devices may include in one embodiment storage 
devices, such as a floppy disk drive 3 50, and input 
devices, such as keyboard 352 and mouse 354. I/O hub 340 
may also be coupled to, for example, hard disk drive 3 56 

10 and compact disc (CD) drive (not shown) . It is to be 

understood that other storage media may also be included in 
computer system 300. 

In an alternate embodiment, the I/O controller 346 may 
be integrated into the I/O hub 340, as may other control 

15 functions. PCI bus 344 may also be coupled to various 

components including, for example, a memory 360 that in one 
embodiment, may be a multilevel, segmented unified memory 
device much as has been described herein. Additional 
devices may be coupled to the I/O expansion bus 342 and to 

20 PCI bus 344. Such devices include an input/output control 
circuit coupled to a parallel port, a serial port, a non- 
volatile memory, and the like. 

Further shown in FIG. 3 is a wireless interface 3 62 
coupled to the PCI bus 344. The wireless interface may be 
25 used in certain embodiments to communicate with remote 

devices. As shown in FIG. 3, wireless interface 3 62 may 
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include a dipole or other antenna 363 (along with other 
components not shown in FIG. 3) . While such a wireless 
interface may vary in different embodiments, in certain 
embodiments the interface may be used to communicate via 
5 data packets with a wireless wide area network (WWAN) , 
wireless local-area network (WLAN) , a BLUETOOTH™ - 
compliant device or system or another wireless access 
point. In various embodiments, wireless interface 362 may 
be coupled to system 300, which may be a notebook personal 

10 computer, via an external add-in card, or an embedded 

device. In other embodiments wireless interface 3 62 may be 
fully integrated into a chipset of system 300. 

Although the description makes reference to specific 
components of the system 3 00, it is contemplated that 

15 numerous modifications and variations of the described and 
illustrated embodiments may be possible. Moreover, while 
FIG. 3 is a block diagram of a particular system (i.e., a 
notebook personal computer) , it is to be understood that 
embodiments of the present invention may be implemented in 

2 0 another wireless device such as a cellular phone, personal 
digital assistant (PDA) or the like. 

In addition, skilled practitioners recognize that 
embodiments may also be realized in software (or in the 
combination of software and hardware) that may be executed 

2 5 on a host system, such as, for example, a computer system, 
a wireless device, or the like. Accordingly, such 
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embodiments may comprise an article in the form of a 
machine -readable storage medium onto which there are 
written instructions, data, etc. that constitute a software 
program that defines at least an aspect of the operation of 
5 the system. The storage medium may include, but is not 
limited to, any type of disk, including floppy disks, 
optical disks, compact disk read-only memories (CD-ROMs) , 
compact disk rewritables (CD-RWs) , and magneto-optical 
disks, and may include semiconductor devices such as read- 

10 only memories (ROMs) , random access memories (RAMs) , 
erasable programmable read-only memories (EPROMs) , 
electrically erasable programmable read-only memories 
(EEPROMs) , flash memories, magnetic or optical cards, or 
any type of media suitable for storing electronic 

15 instructions. Similarly, embodiments may be implemented as 
software modules executed by a programmable control device, 
such as a computer processor or a custom designed state 
machine . 

Accordingly, from the Description above, it should be 
2 0 abundantly clear that embodiments of the subject invention 
constitute a substantial embellishment in spatial audio 
rendering techniques. To wit: an algorithm for spatial 
audio rendering in which filters are applied to simulate 
sound reverberation in a computationally effective manner. 
25 In addition, because that architecture of the spatial audio 
rendering system incorporates a filter bank having 
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parameters that are tunable to a predetermined number of 
reverberation paths, the system facilitates an exercise of 
design discretion in which computational complexity and 
quality of audio reproduction may be balanced. 
5 While the present invention has been described with 

respect to a limited number of embodiments, those skilled 
in the art will appreciate numerous modifications and 
variations therefrom. It is intended that the appended 
claims cover all such modifications and variations as fall 
10 within the true spirit and scope of this present invention. 



23 



